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Abstract 

The detection of an electromagnetic transient which may originate from a binary neutron star 
merger can increase the probability that a given segment of data from the LIGO- Virgo ground- 
based gravitational-wave detector network contains a signal from a binary coalescence. Additional 
information contained in the electromagnetic signal, such as the sky location or distance to the 
source, can help rule out false alarms, and thus lower the necessary threshold for a detection. 
Here, we develop a framework for determining how much sensitivity is added to a gravitational- 
wave search by triggering on an electromagnetic transient. We apply this framework to a variety of 
relevant electromagnetic transients, from short GRBs to signatures of r-process heating to optical 
and radio orphan afterglows. We compute the expected rates of multi- messenger observations in 
the Advanced detector era, and find that searches triggered on short GRBs — with current high- 
energy instruments, such as Fermi — and nucleosynthetic 'kilonovae' — with future optical surveys, 
like LSST — can boost the number of multi-messenger detections by 15% and 40%, respectively, 
for a binary neutron star progenitor model. Short GRB triggers offer precise merger timing, but 
suffer from detection rates decreased by beaming and the high a priori probability that the source 
is outside the LIGO- Virgo sensitive volume. Isotropic kilonovae, on the other hand, could be 
commonly observed within the LIGO- Virgo sensitive volume with an instrument roughly an order 
of magnitude more sensitive than current optical surveys. We propose that the most productive 
strategy for making multi-messenger gravitational-wave observations is using triggers from future 
deep, optical all-sky surveys, with characteristics comparable to LSST, which could make as many 
as ten such coincident observations a year. 
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I. INTRODUCTION 



The first generation of ground-based interferometric gravitational- wave (GW) de- 
tectors was successfully deployed during the last decade. The Laser Interferometer 
Gravitational- Wave Observatory (LIGO 1 ), the Virgo detector 2 , and the GEO-600 detec- 
tor 3 have taken data that were analyzed to search for gravitational- wave signatures of com- 
pact binary mergers, short-lived transients, non-axisymmetric neutron stars, and stochastic 
gravitational- wave backgrounds l|, |2J. No detections have been made so far, as searches have 
resulted in increasingly stringent upper limits on binary merger rates [e.g., [j, Q). Currently, 
the network is undergoing upgrades to the "advanced" configuration, which will increase the 
detector sensitivities by a factor of ~ 10, and thus the expected detection volume by ~ 10 3 

Coalescences of compact-object binaries composed of neutron stars (NS) and black 
holes (BH) are expected to be among the most promising sources of gravitational waves for 
advanced LIGO and Virgo detectors. The expected rates of local GW sources have been 
explored for over two decades [e.g., |7H9|. Binary neutron stars (BNS) are predicted to merge 



at a rate of between 0.01 and 10 coalescences per Mpc 3 of comoving volume per Myr 10 ]. 
The advanced GW detectors will be sensitive to BNS mergers to distances beyond 400 Mpc 
for optimally oriented systems, with a typical range of ~ 200 Mpc. At this sensitivity, even 
conservative predictions indicate that the first detections could be made soon after advanced 
LIGO/Virgo become operational in 2015-2016 flol |. 

Coalescences involving neutron stars are also expected to produce a variety of elec- 
tromagnetic transients. First, and perhaps most notably, NS-NS/NS-BH are promising 



candidates for the progenitors of short gamma-ray bursts (SGRBs) [e.g., |ll|, [l2|. Unlike the 
long-duration majority of GRBs — associated with the final collapse of massive stars [e.g., 



13 1 — SGRBs are believed to be ietted emission resulting from a small fraction of a solar 

r n 

mass rapidly accreting onto a stellar-mass black- hole [e.g., 114J]. In addition to gamma and 
X-ray emission, lower-energy observations have, on rare occasion, been associated with the 



prompt' emission of long GRBs [e.g. [l5|, but so far not with short GRBs [16|]. Short and 



1 http : //www . ligo . caltech . edu 



|http : //www . virgo . inf n . it| 
3 http://www.geo600.org 
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long GRB progenitors have also been associated with late afterglow emission observed from 
the X-ray to the radio [e.g., ll7H22| on timescales between hours and weeks. In addition to 
the prompt and afterglow emissions associated with SGRBs, the neutron-rich tidal-ejecta — 
produced just prior to coalescence — have been proposed as a site of r-process element pro- 



duction 
day 25 



23 



24] . and associated optical emission ('kilonovae') on the timescale of about a 



31] . Throughout this paper, we quote results for the binary neutron star progenitor 
model, using the BNS GW detector horizon distance and inspiral rate; however, our methods 
are equally applicable to NS-BH binaries. 

The prospect of a simultaneous detection of gravitational- wave and electromagnetic 
(EM) signatures from the same event would be particularly exciting. Two channels of 



information from the same source would enable multi-messenger astronomy [e.g.,|32|, probing 
the conversion of gravitational binding energy into electromagnetic radiation. Additionally, 
the detection of a GW signal associated with a SGRB would be a unique and definitive 
determination of SGRBs' cryptic progenitors. The observation of r-process 'kilonovae', from 
a compact binary source confirmed by coincident GW detection, would permit a study of 
the associated densities and pressures through observations of nucleosynthetic products. 
Electromagnetic observations can also yield redshift measurements, which would allow for 
alternative probes of cosmology when combined with distances extracted from GW signals 



33] 



Rather than relying on serendipitous observations of EM and GW signatures from the 
same events, two alternative approaches to increase the detection rate have been proposed: 
(1) following up GW detection candidates with target-of-opportunity searches for EM coun- 
terparts, or (2) triggering searches of archival GW data based on observed EM transients 



34] and references therein). 

Approach (1) has been frequently discussed in the literature [e.g., 



35 



and recently, 



several follow-up searches for electromagnetic signatures associated with possible GW candi- 



dates have been carried out 



(2) 



44 



39 



40] . In the context of searches for binary mergers, approach 



has so far only been applied 4 to searches of GW data based on GRB triggers 



45] as proposed by e.g., Kochanek and Piran (46j and Finn et al. 47] 



4 High-energy neutrinos 42 1 and soft-gamma repeaters / anomalous X-ray pulsars [e.g.. 143] have also been 
used as triggers in searches for GW signals. 
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In this paper, we focus on approach (2) by performing a careful analysis of the im- 
provement in sensitivity when information from an observed electromagnetic transient is 
used to trigger a GW search in archival LIGO/Virgo data. Such information can boost the 
a priori probability that a detectable GW signal exists in the data. EM triggers can fur- 
ther reduce the false alarm rate in GW searches by providing additional constraints on the 
sky location, inclination, or distance to the source. We use a Bayesian framework (section 
HT|) to estimate the amount by which GW search thresholds can be reduced in EM-triggered 
searches relative to optimal, coherent all-sky searches. We go beyond GRB triggers and con- 
sider a set of telescopes and EM surveys (section HTT1) to identify classes of merger-associated 
electromagnetic transients that could increase detection rates. 

We find that, for short GRB triggers (assuming these arise from compact-binary merg- 
ers), the threshold reduction is less than 10% — far more modest than the ~ 50% reduction 

n 

predicted by Kochanek and Piran [461 ] . This difference is due to a more realistic treatment 
of GW detector noise and GRB distribution and beaming angles in the present work, and is 
discussed in Sec. |V] The resulting increase in the total number of GW detections is marginal 
for triggered searches from any EM event, as only a fraction of mergers detectable with LIGO 
and Virgo are expected to have EM counterparts detectable with the concurrent generation 
of EM observatories. The rate of multi-messenger observations, however, can be increased 
by as much as 30-40% for optical, and 15% for high-energy triggers. While the benefits 
of triggering searches on short GRBs will be valuable for determining their progenitors, 
the most advantageous strategy for making multi-messenger observations is using r-process 
kilonovae detections from deep, optical all-sky surveys, such at LSST. 

II. FORMALISM 

Standard, all-sky LIGO- Virgo searches for gravitational waves from coalescing com- 
pact binaries rely on matched filtering against a bank of template gravitational waveforms 
in order to extract signals from the noise j^sj. So far, most searches have not been fully 
coherent between detectors; instead, in coincident searches the best-fitting templates are 



found separately for each detector, and candidates are selected by searching for 
in multiple detectors that are compatible in time and mass parameters [e.g. 



oud events 



49|. This 



approach is suboptimal, since information about the relative phases and amplitudes in dif- 
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ferent detectors is not taken into account, and some candidates may not correspond to a 
self-consistent solution for extrinsic parameters, such as sky location, inclination, and dis- 



tance to the source. Fully coherent searches have been developed [50|, |51f, but these are 
computationally expensive, and have not been regularly used except when an EM trigger is 
used to select the source location on the sky and to limit the time window, which makes a 
coherent search significantly less expensive. 

In this paper, we compare the improvements due to the presence of an electromag- 
netic transient trigger relative to a fully coherent blind all-sky search, in the anticipation 
that all-sky coherent searches will be computationally feasible when advanced detectors are 
operational. Otherwise, if computational costs prevent coherent searches except when the 
sky location is known from the presence of an EM transient, the enhancements due to EM 
triggers will be even greater than we estimate. 

A. Bayes' Rule 

The Bayesian formalism allows us to compare multiple hypotheses based on the given 
evidence and prior beliefs. Consider an observed data set d and a set of competing models 
{Mj|z = 1, 2, . . .} to explain that data, each with continuous parameters Q^. Given a model 
and its parameters, the likelihood of observing the experimental data, L(d\9i, Mi), can be 
predicted. Bayes' rule allows the posterior probability distribution function to be computed 
for a given set of parameters given the assumed model and the observed data: 

p(WM * ) = — rtW) — ' () 

where p(9i\M i ) denotes the a priori probability distribution of the model parameters before 
the data is taken into account. The denominator p(d\Mi) is a constant determined by the 
requirement that posterior probability p(8i\d, Mi) must be normalized to one: 

p{d\Mi)= [ d6iL{d\6 h Mi)p{6i\Mi) , (2) 

where Vi is the parameter space volume in model Mj. This value, known as the 'evidence', 
can be interpreted as the overall probability of generating the observed data set if model Mi 
is correct. 
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Alternately, if we are interested in the posterior probability of a particular model given 
a set of data, p(Mi\d), we can apply Bayes' rule as 

, , n p(d\Mi) p(Mi) , , 

p(Mi\d) = n 1 l !r %) . 3 
p(d) 

Here p(Mi) is the prior for model Mi, and p(d) is another normalization constant defined 
discretely as 

p(d) = J2p(d\M l )p(M l ), (4) 

i 

under the assumption that all alternative models have been enumerated. When comparing 
two alternative models, M\ and M 2 , this normalization cancels, so that the 'odds ratio' 
between the models is just 

Q _ pjM^d) _ p(Mi) pjd\W) 
p(M 2 \d) p(M 2 ) p{d\M 2 y 

where the first term is the ratio of prior probabilities of the models and the second term, 
equal to the ratio of their evidences, is known as the Bayes factor. 

Using this formalism, we can compare a blind GW search with one triggered and 
constrained by an electromagnetic transient. 

B. Gravitational Wave Detection 

In a given set of data, there is some probability that the data contains a (detectable) 
GW signal, p(GW\d); and some probability that there is only noise, p(N\d) = 1 — p{GW\d). 
Specifically, we define p(GW\d) as the probability that there is a GW signal in the data 
ending in a time interval of duration r (ms) from a binary inspiral within a distance 5 
(Mpc) 5 The odds ratio for the data to contain such a signal is 

Q _ p(GW|d) = p(GW) p(d|GW) 
~~ p(N|d) p(N) p(d|N) 



p(GW) 
p(N) 

data collected, primarily determined by the expected merger rate of compact binaries within 



The factor p ^ is the prior probability of the data containing a signal, regardless of the 



5 We use r = 100 ms, corresponding to a typical time window of a coherent Bayesian analysis; and 6 — 1000 
Mpc, which is sufficiently large to ensure that all binary neutron star mergers detectable as GWs fall within 
this volume. The exact values are not important - they are only relevant for making the Monte Carlo 
simulation described in section Til Dl more efficient. 
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the distance 5. The Bayes factor, B = ^jjp^p, on the other hand, is a purely a posteriori 
measure of the confidence in having observed a GW signal: 

= p(d|GW) = f L(9\GW)p(d\6)d6 

p(d|N) p(d\N) ' 1 ' 

Under the assumption of stationary, Gaussian noise, the Bayes factor can be approximated 
as a function of the signal-to-noise ratio (SNR), 

£ocr^ (SNR)2 , (8) 

where rj is the fraction of the prior volume of the parameter space to which the signal's 



parameters can be constrained 6 (see appendix of 52] , where an exponential is missing from 
Eq. (A.5)). Therefore, reducing the prior volume by a particular factor in the region where 
the likelihood is negligible is equivalent to an increase in the evidence integral by the same 
factor. 

Combining Eqs. [6] & El the odds ratio can be expressed as, 

Ooc^.,-e^. (9, 

The confidence in detecting a signal (O) is thus tied to the observed SNR. In practice, 
a 'detection' is defined by a threshold on the odds ratio such that the false-alarm rate is 
reduced to some (small) predetermined value (e.g., < 1CT 2 yr _1 ). 

Although the odds ratio as written in Eq. (Q corresponds to the assumption of sta- 
tionary and Gaussian noise, we know from experience that LIGO and Virgo noise has non- 
stationary, non-Gaussian artifacts ('glitches'). A confident detection may require an SNR 
threshold of ~ 8 in at least two detectors, correspond to a network SNR of 11.3 [e.g., |4|. 
Therefore, in our analysis we choose to use the artificially large odds ratio threshold of 
6 x 10 15 , based on p = 11.3, r) = 1(T 5 (see table]! and p(GW) « = 1(T 7 (using the 

pessimistic merger rate from from Abadie et al. [l0(] 7 ) in order to empirically account for the 
data quality in a conservative manner. The same odds ratio threshold is assumed to hold 
for GW-only candidates and for those with electromagnetic triggers, O gw = O em = 6 xlO 15 . 



6 The fraction of prior volume in which the likelihood is significant, rj, in turn depends on the SNR: higher 
SNR yield better parameter constraints and smaller r\. We neglect this dependence here. 

7 For a binary neutron star merger rate of 0.01 Mpc~ 3 Myr -1 [lOj], with r = 100 ms and (5 = 1000 Mpc, 
p(GW) = 1Q- J . 
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Our analysis is fairly insensitive to the choice of the odds ratio threshold — or, alter- 
natively, the choice of the priors — due to the large uncertainties in other parameters and 
the super-exponential dependence of the odds ratio on the SNR. However, we note that 
a lower odds ratio threshold — as would result from more stationary and gaussian detector 
noise — would make the corresponding SNR threshold more sensitive to constrains on the 
prior and parameter spaces from EM triggers. 



C. Electromagnetic Counterparts to GW Mergers 



Additional evidence for the presence of a gravitational wave signal decreases the re- 
quired signal-to-noise ratio for a given overall detection confidence (odds-ratio). Because 
the SNR is inversely proportional to the source distance for a given type of inspiral event, 
increasing the effective SNR amounts to increasing the detector horizon by the same factor, 
and the accessible detector volume by its cube. Thus, even a small improvement to the 
SNR can have a large effect on the expected detection rate. We can compute the fractional 
increase in sensitivity, i.e., the fraction by which the SNR threshold can be lowered when 
an EM transient trigger has been observed (denoted by the condition \EM), by rearranging 
Eq. ©: 

r In ( o GW 



C 



SNR 



SNR 



EM 



p(GW) 
p(N) ' /gw 



-1 



In O 



P (GW|EM) 
_ P (N|EM) ' ^ EM 



(10) 



If the GW prior is significantly increased by the presence of an electromagnetic counter- 



part, i.e., 



P (GW|EM) P (GW) 



then ( > 1.0 and the required B (i.e., SNR threshold) is 



p(N|EM) " p(N) 

decreased relative to an analysis of the GW-data alone. The presence of an EM signature 
also affects ( through the parameter space term (jj). Due to the higher spatial resolution 



of photonic observations [see, e.g., 



53] , only a portion of the total GW parameter space will 



be compatible with an EM signature. The effects of an EM counterpart on the GW-prior 
and parameter-space terms are examined in the following sections. Note that the division 
of information from an EM transient into two distinct terms, prior and parameter-space, is 
purely conceptual. Constraints present in both terms have the same effect of ruling out false 
alarms and improving the chances of detecting a GW signal. 
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1. GW Prior 



The prior probability of a data-segment containing only noise (i.e. no detectable GW 
signals) can be expressed as p(N) = 1 -p(GW) and p(N|EM) = 1 — p(GW|EM). 8 The prior for 
a detectable GW is given by the expected merger-rate density 1Z, multiplied by the detector 
volume and time-duration of the data segment: 

4"7T 

p(GW) = K ■ — 5 3 ■ t gw . (11) 

In the case of an electromagnetic counterpart, the GW prior must take into account 
the probability that the EM signature was produced by a merging binary (J 7 ), as well as 
the probabilities that the merger took place within the time interval r GW and within the 
distance 5. More precisely, the probability of having a coalescence given a particular type of 
electromagnetic signature is J 7 . In our analysis, we optimistically assume that J 7 is unity; 
however, even assuming T = 0.1 would decrease the multi-messenger detection rates by less 
than 5% for some of the transients. In general, a triggered search yields benefits over a blind 
all-sky search as long as the product of the prior probability for a GW signal and reduced 
parameter volume rj (discussed in the next subsection) is greater given an EM transient 
trigger than without one. This condition yields a threshold for J 7 , which, given a particular 
EM transient, determines whether or not triggering will be beneficial. 

If we denote the temporal localization of the merger based on the electromagnetic 
signature by At, and the electromagnetically accessible volume by V EM , the probability that 
there was a binary merger within the time r GW and distance 5 given the EM transient 
observation is 

P(GW| EMH^.(W (12) 

For a SGRB, r GW /At would be relatively large — as SGRBs are expected to quickly follow 
the merger |38|]. Orphan afterglows (those observed without prompt SGRB signals), on 
the other hand, would be expected on much larger timescales (days to months) — leading 
to a correspondingly smaller prior. In general, higher telescopic precision better constrains 
the parameter space — increasing ( fill C 2j) . For example, if the EM signal has a distance 



determined to be within the GW-horizon, the volume term in Eq. ( TT2|) goes to unity. If, 



It is also possible that, given an EM signal, there could be a chance, unassociated GW-signal detectable 
in a blind search; we do not consider this possibility here. 
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however, there is no distance measurement (i.e., no redshift), higher sensitivity increases the 
accessible volume (V EM ) — and lowers £, because a priori it becomes less likely that these 
transients come from within the sensitive volume of GW detectors. 



2. Parameter Space 

Given a coincident detection of a GW signal and an EM signal from the same source, 
the EM uncertainty in the sky location (and perhaps distance and inclination) will be sig- 
nificantly less than that for the GW signal alone. Therefore, a more restrictive prior can be 
used for GW analyses based on EM triggers. Consequently, the fraction of the prior volume 
for the parameter space of a GW observation, rj, will be enhanced by the EM observation. 
We can decompose rj into distance-inclination (partially degenerate for GW observations 



52| ) and angular terms, 

V = Vd- V<p- (13) 

We can then compare the wide prior volume which must be considered in blind all-sky 
gravitational-wave searches and the restricted prior volume when an electromagnetic coun- 
terpart has been detected. 

Without an EM signature, the angular term may be on the order of ^(GW) = 
ficw/fisky w 1CT 3 , where O nw is the GW detector angular uncertainty — on the order of 



tens of square degrees [52|, l53[; while the distance-inclination fractional parameter space 
might be an order of magnitude less constrained, i.e. ?/d(GW) ~ 0.01. With the presence 
of an electromagnetic counterpart, the angular localization of an EM signature (generally 
arcseconds to arcminutes) is always more accurate than that of a GW antenna: f2 EM < fi GW , 
and thus ?7^(GW|EM) = 1. Similarly, if the EM transient yields a distance determination (i.e., 
a spectroscopic redshift from the signal itself, or an associable host galaxy), ^(GWIEM) = l. 9 
If a distance determination is impossible, the distance parameter space reduces to that of a 
GW-signal alone. These values are summarized in Table fl] 



9 Strictly speaking, even when orphan afterglows enable redshift determination, the inclination of the source 
may be poorly constrained, while observations of SGRB prompt emission will constrain the source incli- 
nation to within the jet opening angle. 
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V4> Vd 1 



GW 


10~ 3 


10~ 2 


10~ 5 


GW EM 


1 


with redshift: 1 
without: 1CT 2 


1.0 
10" 2 



TABLE I: The fraction of the prior volume of the parameter space to which the source can be 
constrained for detections of GWs alone, and GWs with an EM counterpart is decomposed into 
angular and distance-inclination components. An electromagnetic component will always greatly 
improve the angular localization, while the distance will only be better-constrained with redshift 
determination. 

D. Implementation 

We use a Monte Carlo simulation to implement this formalism in determining plausible 
detection properties. For every telescope and transient combination, binary merger events 
are distributed in space according to a fixed merger rate density per unit source time per 



unit comoving volume. We use the default rate of 1 Mpc -3 Myr -1 from Abadie et al. 10 ]. 
The source sky location (9,<p), inclination l, and polarization ip are chosen isotropicaliy. 
Lightcurves (section HLT|) are interpolated to the relevant viewing angle, and detectability is 
determined. The 'observed' properties of the EM transient are then fed into a calculation 
of ( according to Eqs. Ell CEU H21 and Table HI We model redshift as being determined by 
a spectroscopic follow-up if the time-over-threshold (ToT) of the observable EM signal is 
longer than one day and the EM signal undergoes an e-folding in amplitude during this time. 
These requirements are imposed to roughly account for the subset of detected signals which 
warrant followup, and the time needed to obtain it. How accurately the merger-time can be 
determined based on the EM signal depends on how well sampled, and how well modeled 
the light-curve is. As a conservative approximation, we assume the merger-time is only as 
well constrained as the typical time-at-maximum (TaM) of that signal. The SNR of the GW 
signal alone is calculated for a single detector based on the source distance and sky location 
and inclination angles according to the equation 

SNR = 2.0 (1 + zfl G (jpj e(angles), (14) 
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for a source at redshift z and luminosity distance d. The detector horizon distance dn = 445 



Mpc from Abadie et al. 10j is defined as the distance at which the SNR for an optimally 
located and oriented source equals 8.0. The factor of (1 + zf/ & is based on the scaling of the 
waveform amplitude with the (redshifted) chirp mass, and is accurate for low-mass systems 
whose SNR is limited by the bandwidth of the detector rather than the ending frequency of 
the GW signal. The angular dependence 9 (angles) [e.g., |54|, |55|, is given by 

1 1/2 







F 2 (l + cos 2 l) + 4F 2 cos 2 6 



(15) 



with antenna pattern projections, 



1 



F + = - (l + cos 2 9) cos 20 cos 2ip — cos 9 sin 20 sin 2ip, 
F x = - (l + cos 2 9) cos 20 sin 2ip + cos 9 sin 20 cos 2ijj. 



(16) 

included in 



Cosmological distance measures are converted using WMAP-7 parameters [. 
Table IIVI GW detectability is then determined by the threshold SNR > 8 in the absence of 
an EM trigger, and SNR > 8/( when an EM transient is observed. The conversion between 
merger rates in source and observer times is done probabilistically — an event in a year of 
source time has a probability of (1 + z)~ x of being detected in the corresponding observer- 
year. Finally, bootstrapping from the subset of simulated detections, is used to estimate 
parameter variance. 



III. ELECTROMAGNETIC TRANSIENTS 



Gravitational wave sources in the LIGO- Virgo sensitive band (~100 Hz) are dominated 
by neutron-star and stellar-mass black-hole binaries in the final seconds before coalescence. 
Out of the three permutations of source typ es, roughly 10 binary neutron-star systems have 
been observed in our Galaxy [e.g., |57H59|. Four of those systems have sufficiently small 
orbital separations — a < 5i? , or P ^ 0.5 day — that gravitational radiation will merge 
the system within a Hubble time. Exotic formation channels are required to produce such 
systems starting with an initial binary of two massive stars, including two core-collapse 



60- 



Electro- 



supernovae and most likely a phase of common-envelope evolution [e.g. 
magnetic transients associated with GW mergers most likely require the presence of at least 
one NS. If the binary reaches the Roche limit, the NS will be tidally disrupted — ejecting 
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a few hundredths of a solar mass in one or two tidal tails for NS-BH and NS-NS systems, 



respectively [e.g., bj, |27|, M, M, ]Q5\. 
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There is a growing consensus that the expansion of neutron-rich material and r-process 
powered nuclear heating can act as an effective energy reservoir to power fast optical tran- 
sients 



28 



29j . although the precise peak timescale and temperature is dependent on the 
opacity of r-process nuclei, which is currently not well constrained. The most efficient con- 
version of radioactive energy to radiation is provided by those isotopes with a decay timescale 
comparable to the radiative diffusion time through the ejecta. 

The high angular momentum bulk of material from the merger will form a transient 
disk around the existing or newly formed black hole. Cooling is neutrino dominated and pair 



production and/or a relativistic MHD wind can lead to a jetted outflow with T « 100 



while the disk is rapidly consumed on an accretion timescale [69H75j , on the order of a second. 
Because the emitting region must be several powers of ten larger than the compact binary 
that acts as trigger, there is a further physical requirement: the original beamed, relativistic 
outflow would, after expansion, be transformed into bulk kinetic energy. This energy cannot 
be efficiently radiated as gamma rays unless it is re-randomized, which requires relativistic 
shocks. The gamma-rays we receive come from only the material whose motion is directed 
within 1/r of our line of sight — which must lie within the jet angle 9j. At observer times of 
more than about a week, the blast wave has been decelerated to a moderate Lorentz factor, 
irrespective of its initial value. In this 'afterglow' phase, beaming and aberration effects are 
less extreme, emission is observable from a wide range of angles, and is thus sensitive to the 



ejecta geometry 76]. The minimum random Lorentz factor of protons going through the 
decelerating shock is expected to be comparable to the bulk Lorentz factor, while that of the 
electrons may exceed this by a factor of up to the ratio of the proton to the electron mass. 
The energy of the particles can be further boosted by diffusive shock acceleration as particles 
scatter across the shock interface repeatedly, acquiring a power law distribution N(j) oc 7~ p , 
where p ~ 2 — 3. In the presence of turbulent magnetic fields built up behind the shocks, 
the electrons are expected to produce a synchrotron power-law radiation spectrum. 

For an approximately smooth distribution of external matter, the bulk Lorentz factor 



10 Several intriguing scenarios for pre-merger electromagnetic signatures have been proposed [e.g., 
but are not considered here. 
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decreases inversely with time and, as a consequence, the minimum accelerated electron 
random Lorentz factor and the amplified magnetic field also decrease. This implies that 
the spectrum softens in time, leading to late optical and radio afterglow emission. As the 
bulk material decelerates, afterglow emission peaks at progressively lower frequencies. Off- 
axis observers see a rising light-curve reaching a peak when the Lorentz factor drops to 
r(t) ~ l/#obs> followed b y a power-law decrease in luminosity asymptotically approaching 



the on-axis lightcurve 



77H80j. In modeling prompt and afterglow emission, the jet angle can 



be inferred by matching both the observed SGRB rate (assuming BNS progenitors), and the 
observed afterglow luminosities. Based on these req uirements, we use 9j ~ 0.2, consistent 



with jet breaks observed in GRB afterglows [e.g., 



81 



, and references therein]. 



A. Instruments and Surveys 



To explore plausible detection scenarios, we use the parameters of several telescopes 
and surveys across the electromagnetic spectrum. The values used, while consistent with 
each instrument's characteristics, should be taken as representative in an order-of-magnitude 
sense. Some of the instruments addressed have temporary or as-of-yet undecided strategies 
and time-allocations which will decrease their overall time-sky coverage. 



In the optical, we explore both the Palomar Transient Factory [82j (PTF) and Large 
Synoptic Survey Telescope 83] (LSST), which are designed for deep, fast-transient surveys of 
large fractions of the sky. Numerous other optical telescopes are in development or already 



exist (e.g., Pan-STARRS [84]). but have been omitted due to their parametric similarity 
to PTF/LSST, or their focus on longer-cadence observations (e.g., SkyMapper [85|). In 



the radio we examine three surveys — Apertif [86], ASKAP 87|, and LOFAR 88j. Apertif 
is currently taking proposals for survey strategies; out of a range of options, we choose 
(arbitrarily) a very narrow, very deep survey to juxtapose with the shallower, wider ASKAP 
survey. For high-energy observations, we consider the Swift satellite's Burst Alert Telescope 
(BAT) 0] in the X-ray regime (15 - 150 KeV), and the Fermi Gamma-ray Burst Monitor 
(GBM) |90J which extends into gamma-rays (8 KeV - 40 MeV). The parameters used for 
each instrument and survey are presented in Table [TT1 

For detector horizons on the scale of Gpc and larger, cosmological effects become 
important. We assume that the merger rate [lOj is constant in comoving volume until the 
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star-formation-peak at z « 1.5 |9l|, [92| (luminosity distance dh ~ H Gpc), and is negligible 



earlier 



93l | . Detectability is calculated in temporally-redshifted luminosity space, but spectra 



are assumed to be approximately constant between the source and detector frames (i.e. 'K 
corrections' are not considered). 



Project 



Band 



Sensitivity 



FoV (Survey) Cad.(d) Ref. 



X-ray 

Swift (BAT) (15 _ 15Q KeV) 



10" 8 (erg cm -2 s _1 ) 



4,600 
(40,000) 



[89, 94] 



7-Ray 

Fermi (GBM) (g ReV _ 4Q MeV) 



10" 6 (erg cm" 2 s" 1 ) a 



31,200 
(40,000) 



LSST 



Optical 



24.5 



(r: 550 - 700 nm) 5.8 X 10" 30 (erg cm" z s" 1 Hz" 1 ) 



-2 „-l tt -1a 



9.6 

^10,000) 



PTF 



Optical 21.0 7.9 

(r) 1.4 X 10~ 28 (erg cm" 2 s" 1 Hz" 1 ) (8000) 



[82, 95 



Radio 0.1 //Jy 

Apertif ( 1000 _ 175Q MHz ^ 1>Q x 1Q -30 ^ erg cm -2 g -i Hz -i) 



8.0 



ASKAP 



Radio 0.1 mJy 30.0 

(700 - 1800 MHz) 1.0 X 10 -27 (erg cm" 2 s" 1 Hz" 1 ) (20,000) 



Low Radio 1.0 mJy 3,000 

LOFAR ( 1Q _ 20Q MHz ^ 1Q x 1Q -26 ^ erg cm -2 B _i Hz -i) (20,000) 



m 



"This sensitivity, as given in the literature, is specific for the 50 - 300 KeV range. 



TABLE II: Characteristic telescope and survey properties used for analysis of electromagnetic 
transients associated with gravitational- wave progenitors. These values are approximations to the 
true survey designs and strategies, which are, in some cases, yet to be determined. Both the 
instantaneous Field of View (FoV), and the survey FoV — corresponding to the listed cadences 
(cad.) — are given in square-degrees. 



B. Short GRBs 

From a sample of about 60 observed short bursts, roughly 16 have observed redshifts 
determined from spectroscopy of their associated host galaxy. These bursts, compiled by 



Berger 96J are presented in Table IIHI In our simulation, the prompt emission in the BAT 
band is drawn from a log-normal luminosity distribution, constructed to roughly match the 
observed redshift-luminosity distribution of bursts, and the overall BAT event rate, for a 
flux selection cutoff of 10~ 8 erg cm" 2 s" 1 . The luminosity-function parameters are presented 
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in Table IIV} and the distribution is compared with BAT detected luminosities in Fig. [TJ The 
luminosity is extrapolated to the GBM band using the best-fit Band-model of BATSE data 
from Kaneko et al. 97J (see Table II Vp . In both bands, the prompt emission is assumed 



to have a rest-frame duration of one second — approximately the average detected observer- 
frame value for the BAT. The temporal uncertainty connecting SGRBs to GWs, however, is 
taken as six seconds as a more conservative upper limit, and consistent with the assumptions 



made in previous SGRB-triggered GW searches [e.g., [45]. Finally, the emission is assumed 
to be constant within the jet angle, and zero outside. 
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FIG. 1: SGRB luminosity and distance distribution for approximately a year of simulated events 
(grey) and simulated detections, with flux > 10 -8 erg cm -2 s _1 (red). Isotropic equivalent lumi- 
nosity in the BAT band is plotted against luminosity distance and redshift. Over-plotted are the 
SGRBs with known distances from Table ITTT1 (excluding GRB 090426 at z = 2.6). The simulated 
and observed SGRB properties appear consistent, and larger samples of data show populations 
with matching average properties. 
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0.409 


2,240 


0.19 


3.5 x 10 51 


2.1 


< 8.4 x 10 z ' 


061217 


0.827 


5,292 


0.21 


2.7 xlO 51 


2.8 


< 6.7 x 10 28 


070429B 


0.902 


5,896 


0.50 


5.2 x 10 50 


4.8 


< 2.5 x 10 28 


Average 


0.70 


4,707 


0.90 


5.9 x 10 49 


11.7 


2.6 x 10 28 



961 ] . Redshifts have been 



TABLE III: Short GRB and optical afterglow properties from Berger 
converted to luminosity distances, and combined with the T90 — a typical measure of burst duration 
- these were used to convert fluxes to isotropic-equivalent luminosity. The 'average' luminosities 
are medians, and do not include values with only upper limits. All of the X-ray data (corresponding 



to -F 7 in Berger 



96J) correspond to Swift BAT observations, except for GRBs 050709 and 060121. 



C. Afterglows and r-Process Tidal Tails 



While short GRB prompt emission is constrained to the jet-axis, interaction of the 
relativistically jetted material with the circumburst medium produces strong emission at 
much broader angles, and is thus expected to be observable for an off-axis source as an 
'orphan afterglow.' To model these afterglows, we use the synthetic afterglow library of van 
Eerten and MacFadyen [8o|, with a jet energy of 10 48 erg, an ambient density of 1 cm -3 and 
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a jet- angle of 0.2 radians; these values are presented in Table HVl 

Transient Parameters 
Jet Energy (Ej) 10 48 erg 

Circumburst Density (n) 1 cm" 3 

SGRB Jet Half Opening-Angle (0j) 0.2 rad 

Band-Model GRB Spectra a = -1.08, (3 = -2.33, Speak = 262 KeV 

SGRB luminosity function 

PYn ( [Iog(-^/(erg 
erg b-1 ; - CTV ^ eX P ^ 2^- 

Monte Carlo Parameters 
GW Horizon (d H ) 445 Mpc 

GW Data Segment Duration (r GW ) 0.1 s 

Odds Ratio (C EM = C GW ) 6 x 10 15 

Merger Rate (K) 1 Mpc" 3 Myr" 1 

n A = 0.734, n b = 0.0449, n c = 0.222 

Cosmological Parameters" , i , , i 

& H Q = 71.0 km s" 1 Mpc" 1 



p ^ log L ^ _ i. f .i.^L/i.^s .n- ; ,_-^ /' — 10l).7. rr !.(» 



"Larson et al. 

TABLE IV: Afterglow simulation parameters from van Eerten and MacFadyen 



model parameters from 



80], best-fit Band 



971 ]. and additional values used in the current study. 



R-process powered emission from the tidally ejected NS tails is expected to have com- 
parable luminosities to optical afterglows, on shorter timescales for typical viewing angles. 



In our analysis we use the lightcurves generated by Roberts et al. [29|, in their BNS model. 
Figure [2] shows the characteristic light-curves for each type of transient used in our simula- 
tion, along with several observed transients for comparison. It is possible that opacity for 
r-process events is much larger than expected (Kasen, private communication), which would 
decrease their luminosity and drive the main emitted energy to longer wavelengths. Thus, 
our assumptions here may prove optimistic. 
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FIG. 2: Electromagnetic Transients Associated with Compact Binary Mergers. Dots denote peak 

ustration, SGRB 051221 BAT and optical 



2l| are plotted as solid lines, with a single 



luminosities of observed SGRBs and afterglows. For il 
(Gemini-N/GMOS) light curves from Soderberg et al. 
radio (VLA) afterglow observation as an orange dot. Simulated afterglows and r-process tidal tail 
emission (dark green) are shown with dashed lines for different observer angles, as described in the 
legend. Simulated isotropic-equivalent average luminosities and TaMs are plotted with 'x's. 

IV. RESULTS 

The average properties of Monte Carlo events which satisfied their telescope's detection 
criteria are presented in Table IVl Observed and simulated SGRBs and afterglows are plotted 
in Fig. [2j In general, the simulated detection properties are consistent with the distance and 
luminosity distribution of observed SGRBs (Table HIT]) . as reinforced in Fig. [TJ However, note 
that we have assumed a time-at-maximum (TaM) of 6 seconds for short GRBs for consistency 



with the wider time window used in previous GRB-triggered GW searches [4J, |45] to allow 
for a margin of safety. We also find that the average simulated TaM for optical afterglows 
is shorter than for SGRB afterglow observations. 

The fraction by which the SNR detection threshold is lowered for a candidate GW 
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signal associated with the EM transient — the ^-factor — is listed for each transient-telescope 
combination in Table I VI I The mean value is given, along with (ior both the events which 
do, and do not, have redshift determination. It should be emphasized that the mean values 
of £ are averaged over all detected EM events in the simulation — not only those which lead 
to triggered GW-detections, which are most significantly aided by transients with redshift 
determination. The vast majority of events from instruments with large detector horizons 
(e.g., BAT and afterglows with LSST) will occur outside of the LIGO- Virgo sensitive volume, 
which has a BNS horizon distance dn ~ 445 Mpc. Similarly, instruments with detector 
horizons significantly less than the GW horizon (e.g., LOFAR and ASKAP) will not increase 
the number of GW detections regardless of the value of (, because any event within their 
detectable volume will already exceed the blind-detection SNR threshold. 



Mean EM Properties 



Transient 


Telescope 


d L (Mpc) 


Lum (erg s 1 Hz l ) 


Angle (rad) ToT (s) 


TaM (s) 




BAT 


5200 


1.8 x 10 50 (erg s _1 ) 


0.13 






SGRB Prompt 
















GBM 


4800 


6.3 xlO 50 (ergs" 1 ) 


0.13 








LSST 


4100 


1.4 x 10 29 


0.16 


4.0 x 10 4 


1.9 xlO 3 




PTF 


1600 


2.4 x 10 29 


0.13 


9.1 xlO 3 


1.3 xlO 3 


Afterglow 


Apertif 


3200 


2.7 xlO 27 


0.82 


3.6 x 10 6 


2.4 xlO 6 




ASKAP 


130 


3.3 x 10 27 


0.69 


1.7 xlO 6 


1.4 xlO 6 




LOFAR 


8.26 


1.5 x 10 26 


1.0 


2.7 xlO 7 


1.7 xlO 7 




LSST 


460 


3.5 x 10 26 


1.0 


2.4 x 10 5 


1.2 xlO 5 


R-Process 
















PTF 


92 


3.5 x 10 26 


1.0 


2.4 x 10 5 


1.1 xlO 5 



TABLE V: Average properties of simulated electromagnetic transients. Each property is the av- 
erage over detected events for the given telescope. The time-over-threshold (ToT) and time-at- 
maximum (TaM) are calculated in the observer frame, and the luminosity is the average peak, 
isotropic-equivalent. 

The value of ( is increased by a higher prior on the presence of a GW signal, and by 
better constraints on the viable GW parameter space. The best improvement to the prior 
comes from an EM-sensitive volume comparable to (or smaller than) the GW-volume, and 
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from the most precise temporal accuracy. In our model, temporal localizability is determined 
entirely by the time-at-maximum (TaM), while the accuracy in position space is relatively 
constant (see Sec. Ill C2p . unless a redshift (distance) can be determined. The TaM approxi- 
mation is motivated by the uncertainty which would be associated with making only a single 
observation of a given transient, while having accurate models of its source. Positional ac- 
curacy, without redshift, is limited by the precision of GW triangulation (significantly worse 
than any EM localization), and thus constant between transient types. Our requirement for 
redshift determination is based on a time-over-threshold (ToT) longer than a day. Based on 
the ToT values in Table |V] redshift determination is rare in optical afterglows, and never 
occurs for SGRB prompt emission alone. 



c 



Transient Type 


Telescope 


With Redshift 


Without Redshift 


Mean 


SGRB Prompt 


Swift (BAT) 




1.063 


1.063 


Fermi (GBM) 




1.063 


1.063 




LSST 


1.131 


1.000 


1.014 




PTF 


1.135 


1.035 


1.035 


Afterglow 


Apertif 


1.053 


1.000 


1.015 




ASKAP 


1.058 


1.015 


1.025 




LOFAR 


1.034 


1.000 


1.007 


R-Process Tails 


LSST 


1.084 


1.038 


1.056 


PTF 


1.085 


1.038 


1.058 



TABLE VI: Lowered detection threshold factor £. If the criteria for redshift determination are 
satisfied, we assume the distance to the event is fully determined, and the ("-factor is increased. 
The mean value of £ reflects the fraction of simulated EM detections with redshift determination. 

As mentioned earlier, the value of ( alone does not fully determine the benefits of 
triggered searches. It must be considered together with the sensitivity of a given electro- 
magnetic search and the total number of transients expected to be observable with it. For 
example, the highest average value of ( occurs for PTF observations of r-process kilonovae; 
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however, the typical distance for such events is under 100 Mpc — and thus any events within 
that volume are already likely to exceed the blind GW-search threshold. Additionally, the 
expected rate of kilonovae detections with PTF, within that volume, is less than one per 
year. The expected EM transient rates are presented in Table IVHI It is apparent that the 
expected detection rates by LOFAR for afterglows, and PTF for those and kilonovae, are 
too small for much benefit. 

Transient Type Emission Telescope EM Rate (yr -1 ) 

X-ray Swift (BAT) 29.1 ± 0.4 

SGRB Prompt 

Gamma/X-ray Fermi (GBM) 71.2 ± 0.3 



LSST 69 ± 2 

Optical 

PTF 1.1 ± 0.2 



Afterglow Apertif 8 ± 1 

Radio 

ASKAP 5.80 ± 0.01 

Low Radio LOFAR 2.678 x 10~ 3 ± 3 x 10~ 6 

LSST 146.0 ± 0.1 

R-Process Tails Optical 

PTF 0.880 ± 0.002 

TABLE VII: Rates of electromagnetic transient detections produced by our simulations. The 
intrinsic merger-rate density is taken as 10 -6 Mpc -3 yr -1 , in comoving volume and the source's 
rest-frame, from Abadie et al. [lpj ] . 

The expected rates of gravitational-wave detections from EM triggered searches are 
presented in Table IVHH along with the factor increase (gain) relative to the 'blind' (non- 
triggered) detection rate of about 30.8 yr -1 . Observations of kilonovae with a deep, wide 
survey like LSST suggest that the GW detection rate could be boosted by almost 10%, while 
all other transients show gains at a one-percent level or below. These results suggest that 
triggered searches offer small or effectively negligible gains to the rate of gravitational wave 
detections alone. 

On the other hand, EM-triggered GW searches could give a significant relative boost 
to the number of multi-messenger observations. Assuming 10~ 6 BNS mergers per Mpc 3 of 
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GW Detections (yr x ) 



Transient Type 


Telescope 


Triggered 


Gain vs. blind 


SGRB Prompt 


Swift (BAT) 


0.01 ± 0.003 


1.000 ± 0.001 


Fermi (GBM) 


0.06 ± 0.02 


1.002 ± 1 xlO" 3 




LSST 


0.51 ± 0.09 


1.017 ± 0.004 




PTF 


0.00 ± 0.01 


1.000 ± 0.006 


Afterglow 


Apertif 


0.00 ± 0.04 


1.000 ± 0.002 




ASKAP 


0.13 ± 0.07 


1.004 ± 0.003 




LOFAR 


0.00 ± 0.07 


1.000 ± 0.003 


R-Process Tails 


LSST 


2.74 ± 0.02 


1.089 ± 1 xl0~ 3 


PTF 


0.03 ± 0.02 


1.001 ± 0.001 



TABLE VIII: Rate of gravitational-wave detections from an EM-triggered search, and gain factor 
relative to the blind, all-sky GW detection rate of ~ 30.8 detections per year for the chosen merger 
rate. The presence of certain electromagnetic transients increases the prior probability of LIGO 
data containing a detectable signal. Furthermore, information from the EM observation better 
constrains the allowed GW parameter space. Statistical uncertainties from Monte Carlo modeling 
are included. 

comoving volume per year, advanced detectors are expected to make tens of GW detections 
a year. The intrinsic rate of GW and electromagnetic coincident detections, however, is 
much lower — and shows a much more noticeable improvement from triggered searches. 
The rate enhancement for coincident detections is presented in Table IIXI We compare blind 
searches (in which simultaneous detections of EM and GW transients from the same event 
are purely fortuitous) with a combination of blind and EM-triggered searches. Triggers 
made by high-energy observations of SGRBs show improvements of about 15%. Still, with 
both Swift and Fermi, such coincident detections would only be expected once every two 
to ten years. LSST observations of both afterglows and kilonovae show a noticeable rate 
enhancement of ~ 30% and ~ 40%, respectively. In the afterglow case, the boosted rate is 
still just over two coincident events per year for the default BNS merger rate, while that of 
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kilonovae increases from roughly 7 to 10 detections per year. 



GW+EM Coincident Rate 



Transient Type 


Telescope 


Blind 


Blind + Triggered 


Gain 


SGRB Prompt 


Swift (BAT) 


0.0908 ±9x 1(T 4 


0.1046 ± 9xl0~ 4 


1.15 ± 0.02 


Fermi (GBM) 


0.368 ± 0.002 


0.422 ± 0.002 


1.147 ± 0.008 




LSST 


1.80 ± 0.02 


2.31 ± 0.02 


1.28 ± 0.02 




PTF 


0.096 ± 0.005 


1.008 ± 0.005 


1.05 ± 0.07 


Afterglow 


Apertif 


0.0051 ± 6 x 10~ 4 


0.007 ± 0.001 


1.3 ± 0.2 




ASKAP 


4.64 ± 0.05 


4.77 ± 0.05 


1.03 ± 0.02 




LOFAR 


0.0028 ±8x 10~ 4 


0.0028 ± 8 x 10~ 4 


1.0 ± 0.4 


R-Process Tails 


LSST 


7.14 ± 0.02 


9.88 ± 0.03 


1.384 ± 0.005 


PTF 


0.661 ± 0.001 


0.69 ± 0.001 


1.044 ± 0.002 



TABLE IX: The coincident gravitational-wave and electromagnetic detection rates, for both blind- 
search detections made serendipitously, and detections made with a combination of a blind search 
and searches triggered on EM-transient observations. Due to the lower intrinsic probability of 
making blind coincident detections, the rate enhancement from triggered searches is much more 
significant for coincident detections than for GW detections alone. 
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V. DISCUSSION 



In this study, we have outlined a Bayesian framework for evaluating the increased 
observational sensitivity of gravitational-wave detectors in searches triggered by electromag- 
netic transients. To determine the plausible benefits of such triggered searches, we apply 
this framework to a variety of electromagnetic transients associated with binary neutron-star 
mergers, and a series of telescopes and surveys to identify them. These triggered searches 
decrease the required signal-to-noise threshold for a positive detection through the greater 
a priori probability of the presence of a GW signature in the associated LIGO/Virgo data 
and tighter constraints on the parameter space of possible signals. We find that observa- 
tions of r-process kilonovae by a deep, wide-field transient survey like LSST provides the 
maximum benefits — increasing the rate of multimessenger detections by almost 40%. Once 
advanced LIGO/Virgo and LSST are at design specifications, we predict that using such 
triggered gravitational- wave searches could increase the detection rate to about 10 coin- 
cident detections per year. Using optical triggers from both kilonovae and SGRB orphan 
afterglows, along with high-energy triggers from Swift and Fermi could together increase the 
coincident detection rate to about 14 per year, although some of the binary mergers may be 
double-counted as, e.g., Swift/Fermi and LSST triggers. 

This type of multi-messenger astronomy offers tantalizing prospects for probing ultra- 
compact objects, their binary dynamics, and their eventual merger; in addition to possible 
tests of cosmology in the low-redshift universe, and possible insights into the origin of r- 
process nucleosynthetic elements. A SGRB-GW coincident detection might be the most 
exciting prospect, as it would represent a definitive determination of the progenitor to these 
energetic outbursts. In the explicit absence of such a coincidence, the use of triggered 
searches also improves the range at which such a detection could exclude the binary progen- 
itor hypothesis [98 ] . 




While SGRB triggers offer the highest timing accuracy, the uncertainty in their distance 
and the low event rate within the LIGO sensitive volume make them a sub-optimal trigger 
with current X-ray instruments. Keeping in mind that we have assumed a very simple 
luminosity function (Fig. [I]), it suggests that as many as half of SGRBs within a few hundred 
megaparsecs could be undetectable, even if jetted towards the earth. An X-ray telescope 
with the field-of-view of the GBM and a sensitivity an order of magnitude higher than the 
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BAT might significantly boost the rate of detections within the LIGO-volume, and thus the 
benefits to triggered GW searches. 

A factor of ten to one hundred more optical kilonovae than SGRBs could be observed 
from within the GW-detection volume (depending on the assumed SGRB beaming angle). It 
is important to note that the exact peak times and peak temperatures of kilonovae depend on 
the currently unconstrained line opacities of r-process elements. Under the assumptions used 
here, these kilonovae observations could significantly enhance the rate of multimessenger 
detections, despite their low timing precision relative to SGRBs. Surveys like PTF — with a 
sensitivity to about 21st magnitude — aren't able to probe deeply enough to boost the LIGO 
threshold. Based on our findings, the most productive electromagnetic survey to trigger 
gravitational-wave searches would be an optical survey with comparable cadence and sky- 
coverage to PTF, but which is about an order of magnitude more sensitive (i.e. reaching 
about 23rd magnitude). LSST, with a sensitivity of about 24th magnitude, is more than 
sufficient to fully capture events in the LIGO sensitive volume within the LSST survey field of 
view — and optimize the rate of multi-messenger observations. While LSST triggers increase 
the multi- messenger detection rate by almost 40%, the boost to the overall GW detection 
rate is just under 10% because LSST only observes about one-quarter of the sky. 

Searches triggered on EM transients will add the same number of additional detections 
to both the total of GW detections and the number of coincident GW+EM detections. 
However, the fractional gain in the GW detection rate will be smaller than the gain in the 
rate of coincident detections by a factor equal to the sky coverage of the relevant instrument 
over the timescale relevant to the transient being observed. Our results suggest that despite 
the promising prospects for coincident detections, the use of triggered searches would only 
marginally increase the total number of detected gravitational-wave signals, regardless of 
survey strategy. 

The pioneering study of triggered GW-searches, carried out by Kochanek and Piran 
46| . predicted a rate enhancement of about a factor of 3 — significantly larger than that found 
in the current work. The discrepancy between these results is due to the difficulty of extract- 
ing a GW signal from noisy data. An underlying assumption of the Kochanek and Piran 



46| analysis is that the detector noise is both Gaussian and stationary; whereas in practice 



the noise can be correlated and 'glitchy', which then requires a higher detection threshold 
to achieve the same false alarm rate. Additionally, multiple filter templates are required 
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for a search in which the component masses are not known, constituting a so-called 'trials 
factor', which increases the false alarm rate for a fixed SNR threshold. Finally, Kochanek 
and Piran 46[ assume a one-to-one correspondence between SGRBs and GW signals, which 



allows them to ignore data corresponding to times without a SGRB observation. 

In our analysis, we have made specific assumptions regarding the necessary odds ratio 
for confident detection in the presence of glitchy noise, the expected sensitivity of advanced 
LIGO and Virgo, and the binary merger rate. It is worth considering how changes to 
these assumptions might affect our results. As illustrated by the disagreement between 



our results and those of Kochanek and Piran [46], should future detectors achieve better- 
behaved data (improved data quality for a fixed average noise spectrum), the improvement 
from triggered searches would be enhanced. The best improvements from triggered searches 
come from EM telescopes with sensitive volumes comparable to those of the GW detectors. 
Therefore, if the GW-detector noise spectrum were lowered over time, increasing the LIGO- 
Virgo sensitivity, deeper EM surveys would become useful. Perhaps the most uncertain 
parameter is the BNS merger rate. If the binary merger rate were lower than expected, 
the prior GW probability p(GW) would decrease while the conditional probability given an 
EM transient detection, p(GW\EM), would remain the same. Thus, while a lower intrinsic 
merger rate would decrease the overall rate of GW detections, the benefit and importance 
of triggering would be enhanced. 

In this paper, we have attempted to more precisely determine the plausible benefits 
of triggering gravitational-wave searches on electromagnetic transients. At the same time, 
the statistical framework we have formulated for analyzing the expected enhancement from 
using multiple observational channels is completely generalizable to any system of correlated 
observations. Th e sa me techn ique could easily be applied to space-based interferometers, 



such as LISA 
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100l | or NGO [10l| ; or be used as a boost to high-energy par ticle astronomy 



such as with Veritas 1021 ]. or neutrino astronomy — with, e.g., IceCube 103 1 . 
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